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ca  Jo-H*. 

We  qive  a  condition ' which  implies  that  the  trivial  solution,  Uf'O, 

5J.  A 

of  a  class  of  reaction-diffusion  systems  with  homogeneous  Dirichlet 
boundary  conditions,  is  a  global  attractor  for  all  non-negative  solutions. 
In  certain  cases,  this  condition,  which  relates  the  diffusion  matrix  and 
the  domain  to  a  parameter  which  depends  on  the  nonlinear  tern,  signifi¬ 
cantly  improves  similar  conditions  which  can  be  obtained  from  energy 
estimates.  Applications  are  given  to  equations  arising  in  mathematical 


ecology. 


AMS  (MOS)  Subject  Classifications  -  35B30,k.  35B40,  3fiK50,  15K55 

Key  Words  -  Reaction-diffusion  systems,  stability,  boundary  values 
Work  Unit  Number  1  -  Applied  Analysis 


Department  ot  Mathematics,  Heriot-Watt  University,  Riccarton,  Currie, 
Edinburgh  EH14  4AS,  Scotland. 

- / 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAC.,c>- 75-C-00.'4 • 

This  material  is  based  upon  work  supported  by  the  National  Science 
Foundat ion  under  Grant  No.  MCS78-09525. 


SIGNIFICANCE  AND  EXPLANATION 


Reaction-diffusion  systems  are  systems  of  nonlinear  partial  differential 
equations  which  arise  in  various  aspects  of  science  and  engineerinq,  including 
mathematical  ecology.  Such  equations  can  be  used  to  describe  the  evolution  of 
interacting  and  diffusing  species  in  a  bounded  domain,  together  with  the  usual 
assumption  that  there  is  no  movement  across  the  boundary;  however,  when  there 
is  migration  across  the  boundary,  other  boundary  conditions  are  more  appropri¬ 
ate.  It  is  therefore  of  interest  to  study  the  behavior  of  solutions  of  such 
equations  with  a  variety  of  different  boundary  conditions;  that  is,  we  suppose 
that  the  density  of  each  species  at  the  boundary  is  zero. 

The  purpose  of  this  paper  is  to  give  a  condition  which  implies  that  each 
species  in  the  system  moves  to  extinction.  Such  conditions  are  easily  obtained 
from  standard  "energy"  estimates  and  have  been  observed  by  a  number  of  authors. 
We  shall  give  a  condition  which,  for  a  reasonable  class  of  equations,  provides 
significantly  sharper  results.  In  praticular,  under  certain  conditions  we 
obtain  the  global  stability  of  extinction  for  any  domain  and  for  arbitrary 
diffusion  rates.  Applications  are  given  to  equations  which  describe  the 


ecological  interactions  of  predation  and  competition. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


GLOBAL  STABILITY  OK  STATIONARY  SOLUTIONS 
OK  Ri  '.ACT  ION  -  D 1 KKU  S I  ON  SYSTEMS 


Robert  A. 


Gardner* 


I. 

The  purpose  of  this  paper  is  to  give  a  condition  which  implies  that  the 
trivial  solution  U  5  0  of  the  reaction-diffusion  system 

U  -  DAU  +  F  (U ) ,  U  (0, x)  »  U°(x)  , 


U 


3Q*IR 


0  , 


(1) 

(2) 


is  a  global  attractor  for  all  non-negative  solutions  of  (1),  (2).  Our  condition 

relates  the  diffusion  matrix,  D,  and  SI  to  a  third  parameter  which  depends  on 
the  nonlinear  term.  In  certain  cases,  this  provides  a  considerable  improvement 
over  the  conditions  and  results  obtained  from  energy  estimates  by  Conway,  Hoff 
and  Smoller  [3].  These  results  have  applications  to  systems  arising  in  mathe¬ 
matical  ecology. 

We  assume  that  D  =  diag (d, , . . . , d  ),  where  each  d.  is  a  positive  con- 

1  n  l 

stant,  F  (U)  *  (f, (U),...,f  (U)),  A  is  the  m-dimensional  Laplacian,  and  that 
1  n 

U  is  a  bounded  domain  in  IR™  lying  on  one  side  of  its  boundary,  which  we  assume 
to  be  smooth.  The  interactive  term  F  is  assumed  to  be  smooth  and  to  possess 
bounded  invariant  sets  Z  of  the  form  Z  =  n?_jt0,a^l,  for  sufficiently  large 
a^  >  0;  that  is,  if  U°(x)  is  continuous  and  has  values  in  Z  then  U(x,t)r  E 
for  x  e  il  and  all  t  0  for  which  the  solution  U  of  (1),  (2)  exists.  This 
will  be  the  case  if  we  impose  the  condition 


F (U)  •  v  (U)  ^  0  ,  (3a) 

for  all  U  e  3E  and  x  e  (1;  here,  v(U)  is  an  outward  normal  to  31  at  U. 

It  follows  that  there  then  exists  a  unique,  smooth  solution  of  (1),  (2)  defined 
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for  all  t  0  and  with  values  in  E;  see  [21  for  details.  It  will  also  be 


assumed  that  f.  has  the  form 
i 


f.  UJ)  =  U.M.  (U) 
1  11 


(3b) 


.  th 


In  the  context  of  mathematical  ecology,  M.  is  the  local  growth  rate  of  the  i 
species. 

We  first  give  a  condition  which  implies  that  the  plane  {u^  »  0)  is  a 
global  attractor  for  solutions  of  (1),  (2)  with  data  U°(x)  e  E.  Let 

(w)  ~  sup  {f^((.^,...,£^  j,wf  i  • » .  i  ^  3 )  >  ) 

w 

and  define  G. (w)  =  /  g, (s)ds.  Now  let 
1  0  1 


2 

«=  supCG^fwJ/w  :  0  w  a^}  . 

(Note  that  A^  is  bounded  by  an  expression  which  depends  on  a^  and  the 
Lipschitz  constant  of  g^,  which  is  finite?  see  [4,  Lemma].)  If  R  is  the 
radius  of  the  smallest  ball  containing  (1,  our  condition  is  that 


X^d'S  <  1  ,  (5) 

where  y  >  0  depends  only  on  m  if  m  >  1?  (if  m  =  1  y  also  depends  on  g^ 

2  -1 

and  R  d^  ) .  The  proof  employs  a  comparison  technique  introduced  by  Conway  and 
Smoller,  [4],  to  estimate  u^  from  above  by  the  solution  w  of  the  scalar 
equation 

wfc  -  d^Aw  +  gi(w),  w(x,0)  =  u°(x),  w 1 3«x n  “  0  •  (6) 

Next,  it  is  shown  that  condition  (5)  implies  that  zero  is  the  unique  solution  of 
the  steady  state  equation 

0  -  d^Aw  +  gi(w),  w|3n  -  °  ,  (7) 

associated  with  (6),  and  hence,  that  w(x,t)  must  decay  to  zero  as  t  approaches 
infinity.  Thus  if  (5)  holds  for  1  <  i  <  n,  the  origin  must  be  a  global 
attractor.  It  should  be  noted  that  if  A^  £  0,  (which  occurs  in  a  number  of 


examples),  our  condition  is  independent  of  D  and  £1.  In  Section  HI,  we 


apply  the  above  result  to  equations  which  describe  the  ecological  interactions 


of  predation  and  competition 


Acknowledgement .  The  author  would  like  to  thank  M.  Crandall  for  suggesting  that 


(the  theorem  was  originally  proved  with  the 


hypothesis  X.  <  0  instead  of  (5)) 


We  shall  begin  with  a  discussion  of  the  case  m  =  1  and  X.  0.  The  proof 


though  trivial,  provides  a  simple  geometric  interpretation  of  condition  (5)  which 


may  help  to  extend  our  results  to  a  more  general  setting.  The  crucial  step  is 


to  show  that  zero  is  the  unique  non-negative  solution  of  (7).  If  w  is  a  smooth 


non-negative  solution  on  £1  =  (-£,£)  we  multiply  (7)  by  w  and  integrate  from 


zero  to  x  to  obtain 


where  C  =  d.w“'(0)/2  +  G.  (w(0)).  Hence  the  solutions  of  our  boundary  value  prob 


When 


.  <v. 

the  solutions  of  (7)  must  cross  the  z-axis  when  x  =  +1.  Clearly,  there  is  no 
such  smooth  non-negative  solution.  We  shall  prove  a  similar  result  in  several 
space  variables  by  reducing  the  problem  to  one  dimension. 

Theorem.  Suppose  that  conditions  (3)  and  (5)  hold.  Then  lim  u.(x,t)  =  0, 

-  t-K®  1 

uniformly  for  x  f  !J,  where  u^  is  the  i*^  component  of  the  solution  of  (1), 

(2). 


Proof.  We  begin  by  noting  that  (3)  together  with  the  definition  (4)  of  g^ 

imply  that  g.(0)  =  0  and  that  g.(a.)  <  0,  so  that  0  and  a.  are  respec- 
l  ii—  i 

tively  lower  and  upper  solutions  of  (7).  Moreover,  it  is  shown  in  [4]  that  g^ 
is  uniformly  Lipschitz  (and  hence  Holder)  continuous  on  [0,aj,  so  that  there 
exists  a  unique  classical  solution  w  of  (6)  with  values  in  [0,a J  defined 
for  all  t  >_  0. 

The  desired  comparison  u^  w  is  obtained  by  arguing  as  in  [4].  In  par¬ 
ticular,  we  let  z  =  u^  -  w  and  h(x,t,z)  =  f  (U)  -  g^w),  so  that  z  satis¬ 
fies  the  equation 


Zt-dil2  +  h(x,t.Z).  *|t-0  '  0  ■  (8> 

From  (4),  we  see  that  h(x,t,0)  <_  0,  so  that  {z  ^  0}  is  an  invariant  set  for 

(8);  hence  u^  <  w  in  £1  x  . 

We  now  assert  that  zero  is  the  unique  smooth  solution  of 

-d.Aw  =  g .  (w) ,  w|  =  0  ,  (9) 

1  1  du 

with  values  in  [0,aj.  Suppose  that  this  is  not  the  case,  and  let  wQ(x)  ¥  0 
be  such  a  solution.  Let  D  =  {x  e  »m  :  |x|  <  R},  where  R  >  0  is  chosen  such 
that  i!  C  D,  and  define 


w*(x)  = 


wQ(x) ,  x  e  £1 
0  ,  x  e  D\n 


Clearly,  w,  €  H  (D) 


We  claim  that  w^  is  a  weak  lower  solution  of 


-diAw  =  gi(w),  w  3d  =  0  ; 


(9). 


-4- 


that  is,  wt  satisfies  the  inequality 

/  d.Vw  ’VVdx  <  /  q  .  (w  )^dx  , 

D  *  D 

for  all  t  ft^(D)  with  <f  >_  0.  It  clearly  suffices  to  prove  this  inequality 

CD 

for  all  v  t  (D) .  Let  n  denote  the  outward  unit  normal  to  3ft.  Then 

dw^/dn  £  0,  since  ^  0  in  ft  and  wQ  *=  0  on  3ft.  Let  K  *  supp^,  and 

let  H  =  3ft  n  k.  If  we  give  H  the  orientation  induced  by  n,  we  have  that 

/  d Vdw  /dndo  =  /  d.v>dw  /dndo  <  0  , 

3  (ftnK)  10  H  1  0 

since  ‘f  »  0  on  3  k.  Thus 

/  d.W’  -Vw  dx  *  /  d.V^'Vw  dx  +  /  d.W’-Vyrdx 
D  Knft  K\ft 

=  /  d.^dw  /dndo  -  /  d  <Mw  dx 

3 (Knft)  1  0  Knft  1  0 

<  /  d,*>g,  (w  )dx  «  /  d.*g.(w.)dx  , 

"Knft  1  1  0  d  1  1 

since  g^O)  *  0.  Thus  w4  is  an  H*(D)  lower  solution  of  (9)^.  Clearly 

w*  ■  a.  is  an  upper  solution,  so  that  we  may  apply  an  existence  theorem  of 

Dueul  and  Hess,  17),  to  obtain  an  exact  H* (D)  solution  w^  of  (9) D  which 

satisfies  w(  <  <  w*.  Since  is  bounded  and  g^  is  Lipschitz  continuous, 

2+ci 

the  usual  bootstrapping  arguments  imply  that  w^  f  C  (D) ,  so  that  by  an  exis- 
tence  theorem  of  Amann,  [1),  there  exists  a  maximal  C  (D)  solution  w  satis¬ 
fying  <  v  <  w*.  Since  any  rotation  of  a  solution  of  (9)^  is  again  a  solution 
of  (9)^,  it  follows  that  w  is  rotationally  invariant,  so  that  w  satisfies  the 
O.D.E. 

~di^Wrr  *  ”  ^  wr^  “  9^ (w) ,  w(R)  «  0,  w^(0)  “  0  ; 

(the  boundary  condition  at  r  =  0  follows  from  the  smoothness  of  w) .  We  rescale 
this  equation  to  obtain 

-P(wrr  +  wr)  =  gi(w)  *  W<1)  “  °*  w  (0)  *=  0  ,  (10) 


-5- 


Now  multiply  (10)  by  w  and  integrate  from  r  »  0  to  r  =  1 


(We  have  used  the  fact  that  G.(w(l))  «  G.(0)  =  0).  First,  we  shall  suppose  that 


m  >  1.  Since  w(l)  «=  0,  by  Poincare's  inequality  there  exists  a  constant 


C  >0  such  that  C 


(s)ds,  so  that 


By  the  Sobolev  embedding  theorem,  there  exists  a  constant  C  >  0  such  that 


The  last  two  inequalities  together  with  (11)  imply  that 


C,C  (m  -  1 ) / (C  +  1) i  note  that  y  depends 


where 


only  on  m.  Suppose  that  w(0)  =  0.  Then  there  exists  r  with  0  <  r  <  1 


w(r  ).  Multiply  (10)  by  w  and  integrate  from  r  =  0  to 


and  such  that  w 


and  from  r=r  to  r  =  1  to  obtain 


p (w  (l)/s  +  (m 


contradiction.  Thus  we  have  that  w(0)  >  0  and  that 


This  inequality  violates  condition  (5),  yielding  the  desired  contradiction 


If  m  -  1,  we  must  use  a  different  argument.  Let  F(w,r)  *  (z,-p  g.(w)) 


so  that  the  solution  of  (10)  also  satisfies 


when  8  ■  Bq  for  some  p^  <  0.  Since  solutions  to  (12)  are  unique,  solutions 
of  (12)  are  uniformly  bounded  by  the  particular  solution  of  (12)  when  B  *=  8^ 
for  any  P  with  0  P  _>  B^.  Viewing  8  as  parameter,  we  apply  Gronwall's 
inequality  to  obtain  for  any  8^,8,,  between  zero  and  PQ 

(wtO.Bj)  -  w(0,P2))2  +  (ztO.B^  -  z(0,P2))2  1y181  -  Pj2  , 

where  y  depends  only  the  Lipschitz  constant  of  p-^g.  and  on  a^.  Taking 
P2  =  0  and  8^  =  8^,  we  have  that 

Y  1w  (0)  2  (1)  2  . 

Since  m  =  1,  we  clearly  have  that  w(0)  f  0,  The  proof  now  proceeds  as  above. 

We  shall  now  show  that  lim^^wfx.t)  *  0,  pointwise,  where  w(x,t)  is  the 

solution  of  (6).  If  this  is  not  the  case,  there  exists  f  !)  such  that 

lim  sup  w(x  , t)  >  0.  By  [8,  Lemma  3.8],  there  exists  a  constant  K  >  0  such 
t*^®  0 

that  |Vxw(x,t)|  <  K  for  (x,t)  e  ft  x  [  1 ,«° ) .  We  may  therefore  choose  a  sequence 

{t  }  with  lim  t  =  °°,  and  such  that  if  w  (x)  =  w(x,t  ),  then  lim  w  (x)  = 
n  n  n  n  n 

wQ(x)  in  Ca(ft),  0  <  a  <  1,  and  wQ(xQ)  >  0.  We  claim  that  wQix)  is  a 
solution  of  (9).  If  f  CQ(ft),  then 

|  /  [d.A^w  +  <f<3.  (w  )]dx|  =  lim  |  /  [d.A^w  +  <fq.  (w  )]dx| 
ilO  iO  ''in  in 

(1  n  n 

■  lim  I  /  d.^w  (x,t  )dx| 

'  i  x  t  n 
n  ft 

<  lim  d.||wt(-,tn)  ||  2||  *||  2  . 
n  L  L 

Now  multiply  (6)  by  w^,  integrate  over  ft  x  (1,°°),  and  apply  Green's  theorem 
to  obtain 


2 

/  (  w  dxdt  =  lim 

1  ft  1  T-k*> 


/  (d. |Vw|2/2  +  G . (w) ) dx 


t*T 


t=l 


so  that  lim^  ^||wt  (•  ,t)  ||  2  =  0.  Hence  is  a  weak  (and  therefore  strong) 

L 


-7- 


solution  of  (0);  since  g^  is  Lipschitz  continuous  and  wQ  is  Holder  continu¬ 
ous,  we  must  therefore  have  that  is  a  non-trivial,  non-negative  classical 

solution,  yielding  the  desired  contradiction. 

We  obtain  uniform  decay  of  w  as  follows.  If  y(t)  ■  ||  w  ( • ,  t)  ||  *  _,  ,  we 

L 

obtain  the  differential  inequality  y’  ♦  Cy  /  wg.(w)dx  by  multiplying  (6)  by 

ft  1 

w,  integrating  over  ft,  integrating  by  parts,  and  by  finally  applying 

Poincare's  inequality.  Since  wg^(w)  is  a  bounded  function  which  converges 

pointwise  to  zero,  the  Lebesgue  dominated  convergence  theorem  implies  that 

lirn  ^  /  wg . (w) dx  ■  0,  and  thus  the  above  differential  inequality  implies  that 
ft  1 

lim  y(t)  »  0.  Since  w  is  bounded,  we  obtain  l?  decay  for  any  p  >  2, 
t4'*  — 


and  we  obtain  uniform  decay  from  (9,  lemma  3.1). 


□ 


Our  result  can  be  extended  in  several  ways.  The  functions  f^(U)  may  be 
allowed  to  depend  on  x  and  t;  we  need  only  alter  the  definition  (4)  of  g. 


by  taking  the  supremum  of  f^(U,x,t)  over  all  appropriate  values  of  x  and  t 
in  addition  to  the  variables  u^,  j  ^  i. 


We 


may  also  replace  -d^A  by  an  elliptic  operator  of  the  form 


-i  "  “  l  *] 

1  j,k 


k9  /9xi9xj  ' 


where  the  a  *s  are  constants  which  satisfy  the  condition 
J*  * 

i  2 

E,  a  £  C  >  d.|£|  .  We  modify  the  proof  by  first  performing  a  linear  change 
j»*j»xjk  i 

of  variables  x  -*■  x‘  which  transforms  into  -d^A. 

We  can  also  consider  other  boundary  value  problems  when  m  -  1;  (2)  is 

replaced  by  the  mixed  condition 

Ui  +  ^(xjdu^/dn  “0,  x  -  +  t  , 


where  ft  -  and  (3^  0.  For  simplicity,  we  consider  the  case  0 

Suppose  that  the  condition 

,-V2 


||eill-<  max{|gi(w)  (-G.  (w))"1  |  :  0  <_  w  <_  ai> 


(13) 


-8- 


hold:;.  Wo  claim  that  there  arc  no  non-negative,  non-trivial  stationary  solution-:. 

Such  a  solution  must  connect  the  ray  w  «  {T(-l)z  when  x  =  -£  to  the  ray 

w  -  -8  (t)c  when  x  £,  where  w  *  z.  Condition  (13)  ensures  that  the  former 
i  x 

ray  lies  above  those  trajectories  which  coincide  with  the  curve  OAK  in  Figure  1 
and  that  the  latter  ray  lies  below  those  trajectories  ..hich  coincide  with  the 

-1  / 

curve  OCD  in  Figure  1,  since  the  slope  of  OAK  and  OCD  is  _+  g^ (w)  (-G .  (w) ) 
when  G^  (w)  /  0,  and  the  appropriate  right  or  left  hand  limit  of  this  expression 
when  G^(w)  «=  0.  It  should  be  noted  that  condition  (13)  is  independent  of  D 
and  ft . 

The  above  remarks  indicate  that  our  theorem  should  be  true  in  somewhat 

greater  generality.  For  example,  it  is  reasonable  to  consider  uniformly  elliptic 

variable  coefficient  operators  of  the  form  L.  =  -E  .  ,  ,  (x,  t)  3"”/3x.  3x  , ,  or 

i  3 # k  3,k  l  ] 

mixed  boundary  conditions  of  the  form  (13)  in  several  space  variables.  However, 
the  proof  given  above  doesn't  seem  to  be  adaptable  to  such  situations.  It  seems 
likely  that  a  proof  might  be  found  which  avoids  using  the  extension  procedure  and 
the  rotational  symmetry  of  the  Laplacian. 

III. 

In  this  section,  we  shall  consider  a  few  examples  arising  in  mathematical 
ecology.  Most  attention  has  been  given  to  the  Cauchy  problem  and  to  the  initial¬ 
boundary  value  problem  with  homogeneous  Neumann,  (or  "no  flux")  conditions,  [4], 

{51,  [6J.  However,  Dirichlet  or  mixed  boundary  conditions  may  be  more  appropri¬ 
ate  if  there  is  migration  across  the  boundary.  It  is  worth  noting  that  our 
result  is  not  valid  for  the  pure  Neumann  or  Cauchy  problem,  since  in  these 
cases,  there  may  exist  non-trivial  constant  rest  states. 

EXAMPLE  1.  PREDATION. 

Let  v  be  the  population  density  of  a  predator  species,  and  let  u  be  the 
density  of  its  prey.  Let  M  and  N  be  the  growth  rates  of  u  and  v  respec¬ 
tively;  then  u  and  v  satisfy  the  system 
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The  interaction  is  characterized  by  making  certain  assumptions  about  the  algebraic 
signs  of  M  and  N  and  their  partial  derivatives;  for  a  complete  discussion,  see 
(4],  [51,  [61.  We  shall  assume  that  (i),  M  0  and  N  <_  0  if  u  and  v  are 
both  near  zero,  (ii),  M  <_  0  if  u  is  large  and  N  <_  0  if  v  is  large,  and 
(iii),  Mv  £  0  and  _>  0.  Condition  (i)  requires  that  there  be  a  critical 
population  density  below  which  each  species  goes  extinct  and  (ii)  is  a  resource 
limitation  condition.  The  predator  prey  relationship  is  characterized  by  (iii) . 
One  such  example  is  the  Rosensweig-MacArthur  equations,  where  we  take 

M  =  u(-v  +  6 (u  -  a) (b  -  u)) 

N  =  v(-v  +  cu  -  d)  ; 

a,b,c,d,  and  6  are  positive  constants.  The  phase  diagram  of  the  vector  field 
(uM,vN)  is  given  in  Figure  2. 


Figure  2. 


•J**v 


anm.  <*%?**% 


Wo  have  that  q^(w)  -  iW(w  -  a)  (b  -  w).  A  simple  computation  dhows  that 
2 

Xj  *  M  (a  ♦  b)  /9  -  ab/2] ,  so  that  X^  <_  0  if  a  b/2.  In  this  ease,  u  must 

t ond  to  *e to  independently  of  P,«,  and  fi.  Since  N  <  0  if  n  is  sufficiently 
small,  wo  must  have  that  v  tends  to  r.ero  also,  (even  though  X ^  is  always 

m 

positive  here). 

EXAMPLE  2.  COMPETITION. 

Wo  now  let  u  and  v  be  the  population  densities  of  two  competing  species. 

Those  variables  will  again  satisfy  a  system  of  equations  of  the  form  (14).  How 

ever,  wo  now  replace  hypothesis  (iii)  with  (iii)’,  M  '  0  and  N  *.  0.  (We 

v  —  u  — 

still  assume  that  (i)  and  (ii)  hold).  The  phase  diagram  of  an  example  of  such 
a  system  is  indicated  in  Figure  3. 


For  definiteness,  let 

M  ■  u(-v  +  ^(u  -  AjMbj  -  u)) 

N  -  v(-u  5^(v  -  a^)(b,  -  v) )  . 

As  in  Example  1,  we  have  that  the  oriqin  is  a  global  attractor,  provided  that 
aA  *  b^/2,  i  -  1,2. 
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